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The paper is dedicated to Mihail Samoilovich Livshits, who was first to consider 
a nonselfadjoit operator as a part of an extended sefadjoint scattering system. 

Abstract. We consider a one-body spin-less electron spectral problem for a 
resonance scattering system constructed of a quantum well weakly connected 
to a noncompact exterior reservoir, where the electron is free. The simplest 
kind of the resonance scattering system is a quantum network, with the reser- 
voir composed of few disjoint cylindrical quantum wires, and the Schrodinger 
equation on the network, with the real bounded potential on the wells and 
constant potential on the wires. We propose a Dirichlet-to-Neumann - based 
analysis to reveal the resonance nature of conductance across the star-shaped 
element of the network (a junction), derive an approximate formula for the 
scattering matrix of the junction, construct a fitted zero-range solvable model 
of the junction and interpret a phenomenological parameter arising in Datta- 
Das Sarma boundary condition, see |13) . for T-junctions. We also propose 
using of the fitted zero-range solvable model as the first step in a modified 
analytic perturbation procedure of calculation of the corresponding scattering 
matrix. 
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1. Introduction 

A typical quantum resonance scattering system is composed of a compact inner 
region surrounded by barriers and an exterior reservoir, where the quantum dy- 
namics is free. These components are weakly connected due to tunneling across 
the barriers or via a narrow connecting channels. Non-compact quantum networks 
(QN) are typical resonance scattering systems. Manufacturing of QN with pre- 
scribed transport properties is now a most challenging problem of computational 
nano-electronics. While physical laws defining transport properties of the QN are 
mostly represented in form of partial differential equations, the direct computing 
can't help optimization of design of the QN, because it requires expensive and re- 
source consuming scanning over the space of physical and geometrical parameters 
of the network. The domain of scanning could be essentially reduced in the case 
when there exist an approximate explicit formula connecting directly the trans- 
port characteristics with the parameters defining the geometry and the physical 
properties of the network. 

We derive an explicit approximate formula for the scattering matrix of a simplest 
QN- a junction,- consisting of a vertex domain - a quantum well - connected to 
the outer reservoir decomposed geometrically into a sum of cylindrical leads. The 
corresponding model Hamiltonian is obtained based on Glazman splitting, [3], 
C — > La © Ia of the original Hamiltonian, depending on the Fermi level A, into the 
sum of two operators with complementary branches of the continuous spectra. The 
model proves to be fitted because the corresponding model Dirichlet-to-Neumann 
map ( DN-map, see [73] ) serves a rational approximation of the DN - map of the 
non-trivial component La of the split system. 

In an important alternative class of the resonance scattering systems, repre- 
sented by the Helmholtz resonator, the reservoir can't be decomposed into simple 
components similar to the cylindrical leads, but the finite leads connecting the 
compact subsystem - the resonator - with the reservoir, admit a similar decom- 
position. Then again, a fitted solvable model can be constructed, see [35], based 
on the splitting of the spectral channels in the leads. The model obtained can 
serve again as a first step - a jump start - of the corresponding analytic perturba- 
tion procedure. We postpone the discussion of the Helmholtz resonator and other 
systems with nontrivial reservoirs to oncoming publications. 

Main difficulty of analysis of resonance scattering systems of both above kinds 
on the networks is defined by presence of the eigenvalues of the isolated compact 
subsystem, embedded into the continuous spectrum of the reservoir, separated 
from the compact subsystem. Indeed, for a selfadjoint operator Aq in the Hilbert 
space E, with discrete spectrum, and small eV ,the self-adjoint perturbation, | 
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II < £, defines, for each simple isolated eigenvalue of Aq 
2e < mint^s|As - Xt \ = Ps 

a branch of eigenvalues of the perturbed operator :— Aq + eV represented 
in form of a geometrically convergent series 

A^ = A° + eA°(l)+s2A0(2) + £3A°(3) + ..., 

and the corresponding branch of eigenfunctions, see [28] . 

This standard analytic perturbation approach is not applicable, generally, to opera- 
tors with eigenvalues embedded into the continuous spectrum, in particular to non- 
compact QN where the "spacing" p is zero. Development of radio-location during 
WWII required analysis of scattering problems on the networks of electromagnetic 
wave guides, in particular on junctions. The scattering on the junction is a typical 
perturbation problem for embedded eigenvalues. The perturbation of the problem 
causes the transformation of real eigenvalues on the vertex domain of the junction 
into complex resonances. This problem can't be solved by methods of the standard 
spectral theory of selfadjoint operators. In the paper [40] M.S. Livshits proposed an 
elegant approach to the problem of transmission of electro-magnetic signals across 
the junction, taking into account only oscillatory electro-magnetic modes in the 
wave-guides and neglecting the "evanescent"- exponentially decreasing modes. He 
reduced the calculation of the scattering matrix to calculation of the correspond- 
ing characteristic function and found a real wave conductance for the oscillatory 
modes and pure imaginary wave conductance for evanescent modes. The discovery, 
based on M.S. Livshits ideas, of the connection between the scattering matrix and 
the characteristic function of the corresponding non-self-adjoint operator, see [T], 
was an extraordinary achievement and became eventually a source/basement of a 
series of important results in the theory of the functional models of the dissipative 
operators, see [39l [50l [49] . The approach to the perturbation theory developed in 
these papers permitted to understand the spectral nature of the resonances, but it 
does not help practical problem of optimization of design of quantum or electro- 
magnetic networks with prescribed transport properties. One more detail in the 
above paper [40] was important in this respect. M.S. Livshits completely disre- 
garded the evanescent waves in the wave-guides, which was usual in the papers of 
engineers and physicists on the electro-magnetic wave-guides. But he emphasized 
in [ID] importance of accurate elimination of the evanescent waves. This was not 
done till now. We see now, that the absence of analysis of the evanescent waves 
prevented M.S. Livshits from establishing, at that time, an effective connection 
between the geometry of the junction and the transmisson/reflection coefficients, 
see also our comments below, section 3. 

In the case of operators with continuous or dense discrete spectrum one can sub- 
stitute the Hamiltonian Aq of an unperturbed system by a fitted solvable model 
A"^ , and then develop an analytic perturbation procedure between the perturbed 
Hamiltonian A^ and the model A^. This two-steps idea Aq A"^ ^ A^ of the 
modified analytic perturbation procedure was suggested, in implicit form, by H. 
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Poincare for relevant problems of celestial mechanics, see [62], and formulated in 
an explicit form in 1972 by I Prigogine. In 1972 I. Prigogine, [53], declared im- 
portance of the search of a general practical algorithm for the two-step analytic 
perturbation procedure 

^0 — > — ' 

implementing the above Poincare idea. Prigogine attempted to find an Intermedi- 
ate operator as a function of the unperturbed operator = $(Ao), and he 
wanted to have the above two step analytic perturbation procedure on the whole 
Hilbert space. The search of the corresponding "intermediate operator" contin- 
ued for almost 20 years, but did not give any results. Finally Prigogine declared 
that the intermediate operator with the expected properties does not exist and 
can't be constructed. 

We guess now, that I.Prigogine's suggestion based on the Intermediate Operator 
was very close to success. The idea of Prigogine was commonly used by physi- 
cists in form of effective Hamiltonian of a complex quantum systems, and, after 
essential modification, in for "geometrical intergration" in dynamical problems 
of classical mechanics. 

In our recent papers [H] |43l [27], see also an extended list of references below, 
we suggested a method of accurate elimination of the evanescent waves based on 
the idea of the intermediate Hamiltonian. Our method also permits to accurately 
eliminate the evanescent waves in the case studied by M.S. Livshits. In this paper 
we provide, following the quoted papers, a review of the corresponding modified 
approach to the analytic perturbation procedure and describe, based on [59], an 
algorithm of construction of the solvable model and the procedure of fitting. We 
developed the corresponding general approach to the spectral problems with em- 
bedded eigenvalues in the series of papers [Ml HZl HZl [5HI 1331 [SHI HI] ■ Contrary 
to the original Prigogine's idea, we do a couple of changes: 

1. We search for the Intermediate Operator - the "jump start", see [5S]- A^ , on the 
first step of the above mentioned two-step procedure, not among functions $(Ao) of 
Aq, as I. Prigogine suggested, but among weak (finite-dimensional) perturbations 
of the non-perturbed Hamiltonian Aq, which is close to the method suggested by 
Livshits. 

2. We restricted our analysis to the part of the unperturbed operator on a spectral 
subspace which corresponds to some "essential spectral interval" , contrary to I. 
Prigogine who attempted to find a global Intermediate Operator on the whole 
space. Thus we develop our modified analytic perturbation procedure locally. A 
similar requirement of locality is applied in [38j on the space of initial data of the 
structure-preserving model. 

We begin with two classical examples of the resonance scattering systems, to reveal 
typical difficulties arising from the very beginning when considering perturbations 
of systems with eigenvalues embedded into the continuous spectrum, and discuss 
nearest prospects of the perturbative analysis of these systems. 

Example 1: Helmholtz Resonator. 
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Helmholtz resonator was probably the first resonance scattering system discussed 
mathematically, see [SS]. It is composed of the typical details: the inner domain 
f2i„t, the shell flsheii, and the reservoir flout separated by the shall from flint- 
Consider the Helmholtz equation — Am = Am in a domain e i?3 represented as a 
sum of two disjoint parts fl = flint U flout and a shell flsheii with a small opening. 
Kirchhoff suggested to substitute the problem by the model where the opening is 
pointwise, so that there exist only one common point a e flint U flout U flsheii ■ 
Kirchhoff suggested an Ansatz for the Green- function G\{x,y) of the Helmholtz 
equation in fl with Neumann boundary condition 

dG 

-AGxix, y) - AGa(x, y) = Six - y) 



0,x,ye fl. 



an 



in the form of a linear combination of the Green functions G'™*(x,?/), G'^*{x,y) 
of the inner and the outer problems: 



-AGr'°"*(x,y) 



AG;^ K^^ V) 



5{x 



y) — 



0. 




Figure 1. Hemlholtz Resonator with a point-wise opening at 
the point a and the enlarged detail of the resonator with a narrow 
short channel, 5 « H « A~^/^ 



G\{x,y) - 



Gr*(x,y) -I- A°^'Gl^\x,a), x,y^flout 
A™* Gf{x, a), if ye flout. X e fl^u 



with undefined constants -the Kirchhoff coefficients yl™*. This Kirchhoff 

Ansatz satisfies the equation and the Neumann boundary conditions everywhere 
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on dQ, except the point a, where the Ansatz is singular. The problem of choice of 
the KirchhofF constant and other interesting problems concerning the resonator, 
see [6 5) remained open until recent time, see in this connection the preprint [lOj . 

Example 2: Zero-range potential 

In 1933 E. Fermi, |18j . considered the problem of scattering of neutrons n by 
the nucleon S of Sulphur. He suggested to choose for this problem the model 
Hamiltonian in the form of Laplacian in -^2(^3) defined on smooth functions u € 
L2{R3,) with a singularity at the origin 

^ ' 47r|a;| 

and a special boundary condition imposed on the asymptotic boundary values 

A" = 7 = 7. 

The Laplacian is symmetric and even self-adjoint with this boundary condition, 
and admits explicit calculation of the eigenfunctions: this model is "solvable". 
Fermi suggested to "fit" this model choosing 7 — —Airp'^'^, if — Pq is a small negative 
eigenvalue in the system n, S. 

The model can be extended to the case when 7 > 0, and fit to the purely imaginary 
resonance pq — i^. The resulting mysterious "zero-range potential" suggested by 
Fermi was interpreted by F. Berezin and L. Faddeev [7] in terms of von Neumann 
Operator Extensions Theory, [47j. Later this "zero-range potential" was used in 
numerous physical and mathematical papers and books, see for instance [Ij. 
In both above examples the reservoirs are either a large exterior domain, or the 
whole space with single point a; = removed. The first example was also treated 
by the operator extension methods in [17| . producing a zero-range solvable model 
of the resonator immersed into 3D space. The role of the unperturbed operator in 
|17j played an orthogonal sum of the Neumann Laplacian Lint in L2{^int) and Lout 
in L2{^out)- The basic difficulty of the original perturbation problem, with a thin 
short channel, is caused by presence of the eigenvalues of Lint embedded into the 
continuous spectrum of Lout- The standard self-adjoint spectral theory is generally 
unable to treat the problem of embedded eigenvalues, by observing transformation 
of them into the corresponding complex resonances. The elegant Lax-Phillips ap- 
proach to resonance scattering problems reveals the spectral nature or resonances, 
see [39], but does not help to calculate them. In [17] the resonances can be easily 
calculated via solving an algebraic equation, but yet the fitting of the suggested 
zero-range model remained a problem. In jlOj an approach to the problem of fitting 
of the model is suggested based on an explicit formula connecting the "full" scat- 
tering matrix of the Helmholtz Resonator with the Neumann-to-Dirichlet map, see 
[54) . and a subsequent rational approximation of the Neumann-to-Dirichlet map 
for the inner domain of the resonator (the cavity ^int)- Fortunately the problem of 
search of the resonances, in the case of small opening, becomes finite-dimensional 
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after replacement of the Neumann-to-Dirichlet map of the cavity by the corre- 
sponding rational approximation, see more details in |10j . 

In the second example just a self-adjoint operator — is suggested, with only 
parameter 7, which can be interpreted in spectral terms. This operator plays a 
role of an effective Hamiltonian of the original scattering problem for the neutrons 
and the nuclei, see [18]. Yet again, the substitution of the original perturbed {full) 
Hamiltonian by the effective solvable Hamiltonian — A-y requires fitting of the 
model, at least on some essential interval of energy. 

Note that the role of the effective Hamiltonian is played, in the second example, by 
a self-adjoint extension of the unperturbed Hamiltonian —A. Numerous effective 
Hamiltonians in quantum mechanics are constructed as zero-range solvable models 
of quantum systems see for instance (53i and our recent papers quoted above. 
In this review we represent some results of our recent papers quoted above (see the 
text preceding the Example 1) where the effective Hamiltonians are constructed as 
zero-range solvable models. To make the text easily readable, we omit some proofs 
and most of technical details, which can be found in the original publications. But 
we pay additional attention to the interconnections of our constructions previously 
spread in different publications. 

2. Scattering on Quantum Networks 
and Junctions via DN-map. 

The basic idea of analysis of partial differential equations on quantum networks 
is that the corresponding Schrodinger problem can be divided in two parts: a 
Schrodinger equation on the region surrounded by barriers ( e.g. a quantum well 
) and one on the reservoir the two being weakly coupled by tunneling, see for 
instance |63j or by a thin channel. It is noticed in [63j that this decomposition 
"corresponds to the schematization of the transport process as a coherent process 
on the quantum well, fed by the exterior reservoirs - quantum wires" . On the 
reservoirs, assumed to be homogeneous and neutral, the electronelectron interac- 
tion is neglected, and the single electron is free. But the resonance properties of 
the quantum well and the tunneling on the contacts define the transport proper- 
ties of the whole network. It is a common belief that thin quantum network can 
be modeled by a Id graph, see |51| . with either Kirchhoff boundary conditions, or 
just non-specified selfadjoint boundary conditions at the vertices. There exist an 
extended bibliography concerning one-dimensional models of quantum networks, 
see for instance [121 |36l [37l [25l [72] . Notice that even sharp resonance effects on 2d 
wave-guides and networks were studied theoretically mainly by numerical methods, 
see for instance [531 US] ■ 

Quantum network which is being studied in this paper, is composed of straight 
leads (quantum wires) width S, some of them semi-infinite, and vertex domains 
ris ( quantum wells), see Fig[2l An important basic detail of the quantum network 
is a junction, see Fig. [3] The junction is a non-compact quantum network com- 
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Figure 2. Quantum Network: a detail 



posed of a quantum well and few semi-infinite quantum wires, of constant width, 
attached to it. The junction is usually called thin, if the diameter of the quantum 
well riint strongly dominates the width S of the wires uj attached to it: S << diam 
f2i„t. Calculation of the scattering matrix of a junction is a challenging computa- 
tional problem, yet accessible for standard commercial programs, see the discussion 
below. Physicists have certain preferences about the boundary conditions at the 
vertices, see the discussion below. Example 3. 

Example 3: Thin symmetric T-junction. For thin symmetric T-junction, with the 
"bar" orthogonal to the "leg" , a reasonably simple explicit formula for the scatter- 
ing matrix was suggested in jl3j based on reduction of the 2D scattering problem 
on the junction to the corresponding ID scattering problem on the correspond- 
ing quantum graph, see Fig. [5] The boundary conditions for the model T-junction 
suggested in [13] , is presented in terms of limit values of the wave- function on the 
ID wires {V'i}i=i ■= "0 and the values of the corresponding outward derivative ( 
boundary currents) ■= '0' a-t the node: 



(1) 



or in the form 







(2) 



with the projection 




(3) 
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Figure 4. Simplest symmetric T-junction with a square vertex domain 



The scattering matrix of such a junction is constant Sp — I — 2P^, with the 
phenomenological parameter (3 responsible for connection between the bar and 
the leg. This formula was intensely used, see for instance [551 [7T], despite unclear 
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Figure 5. Model symmetric T-junction 

meaning of the parameter (3. See further discussion of transmission across the 
junction in [HI [16l [2l| and find more references therein. 

In this paper, based on the resonance conception of conductance on the junction, 
we suggest a semi-analytic procedure of calculation of the scattering matrix and a 
method of reduction of a general thin junction to a quantum graph. Moreover, we 
suggest a solvable model of a thin junction and reveal the meaning of the projection 
P/3. In this paper we do not take into account the spin-orbital interaction, just by 
disregarding the spin of the electron. 

We consider a junction constructed of a few straight leads w™, U^^^w™ — 
u!, width 6, attached orthogonally to the flat pieces of the piecewise-smooth 
boundary of the vertex domain Qint, ^ = ^int U w. On smooth functions u G 
VF2 (^) satisfying the homogeneous Neumann boundary condition, we define the 
Schrodinger operator 

- Au + Vu =: C 

with the potential V equal to the constant Vs on the leads and equal to a real 
bounded piecewise-continuous function on Qi„t supplied with Dirichlet boundary 
condition. The operator C is essentially self-adjoint, and it can be considered as 
a perturbation of the corresponding operator Cq defined by the same differential 
expression and an additional Dirichlet boundary condition on U^^iTm =: F: 

£ > © Lint = £-0 ■ 

The spectrum a{Lint) of Lint is discrete, and the spectrum of L'^ is abso- 
lutely continuous, consists of a countable set of branches a'^ — ^m=i ^hLi '^T 
corresponding to the parts ^™ of l'^ 

with the homogeneous Dirichlet boundary condition u\-^ = at the bottom sections 
X™ L = 0. The operators LJ" on the wires are obtained from L^ via separation of 
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variables, with the basis of cross-section eigenfunctions {e™} = |y^sin ^| , I = 
1, 2, . . . , TO = 1, 2, . . . M. Here the local transversal coordinate on lo"^ is denoted 
by y. The eigenfunctions of L'^ are scattered waves on each lead w™: 

represented as linear combinations of oscillating exponential modes 

with the reflection coefflcients Si — 1. The perturbed operator £ is obtained from 
Co by replacement of the homogeneous Dirichlet boundary condition on the bottom 
sections F by the smooth matching condition. The corresponding scattered waves 
are obtained via matching on F a solution of the Schrodinger equation on the 
vertex domain with the scattering Ansatz (see for instance [45j|46]): "01" (2:) = 

composed, for given A, of the above oscillating modes x!J: the open channels, with 
the thresholds below A, A™ < A, and the exponentially decreasing ("evanescent") 
modes in the closed channels 

= er(y)e-^^-:= e-^-'°%r(y), A < A., (6) 

associated with the thresholds A™ = tt^ 5^^ of the closed channels in the leads 
- see [43] for details. 

The quantum wells and the quantum wires are usually manufactured as a certain 
relief of the surface of the semiconductor. We assume in this paper that the scaled 
Fermi level A — 2m* Epfi^^ of the semiconductor is situated in the middle of the 
first spectral band Ai ^[Vs + ^,Vs + 4f^] of the wire, A = + Then the 
first spectral band plays the role of the conductivity band and the junction has 
metallic properties. At low temperature T, the scattering processes are observed 
only on the essential spectral interval 

= [A - 2m*KTtr^, A - 2m*KTh'^] C Ai. (7) 

If the electron's density is low, the main contribution to the scattering picture 
is defined by one-body processes. In this paper we focus on one-body scattering 
on the essential spectral interval. We disregard the spin-orbital interaction and 
neglect all effects connected with electrons spin. Hence we study the scattering on 
the first spectral band Ai = [vr^ (5^^, 47r^ S^'^] of the open channel, and represent 
the cross-section space L2{T) —: E as an orthogonal sum of the entrance subspaces 
E± of the open and closed spectral channels respectively: 

M M 00 

E+ = y eT, = V V Pe^ =■■ P± ■ (8) 

m— 1 m— 1 1—2 

The infinite linear system for the coefficients of the scattering Ansatz, obtained 
from the matching conditions, can be solved, if the Green functions Gp G,;„t 
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of the Schrodinger operators — Lint in L2{^int), with Dh'ichlct boundary 
conditions is constructed. The operator Lp is defined on W|-functions in rii„t, 
with the Meixner conditions at the inner corner points: 

Lintu — —Am + Vu ~ Xu, u\^^^ = 0. (9) 

The Green function is found from the equation: 

L^G^ = -AG^ + VG^ = XG^ + S{x-y), G^\g^^^^ =0. (10) 

Hereafter we denote by the spectrum of ip . According to the general theory 
of second-order elliptic equations, the solution u of the boundary problem 

— Au + Vu — Xu, 7i|p = Ur, u\gQ^\Y ^ ^' ^^^^ 

is represented by the Poisson map 

u{x) = / Vr{x,j,X)ur{j) d'f , 



involving the kernel Vint{x,j) — —dGp{x,'-f)/dnj. The corresponding boundary 
current on T is calculated as 



du 
dn 



d^G^{x,-f,X) 

-i — Mr (7) dT -. VMr[X)ur 

orixOn^ 



This formal integral operator is restriction onto P of the Dirichlet-to-Neumann 
map, see [TSj |20, ,21,. For the sake of brevity we call it here "relative DN-map". 
The relative DN- map is also a Nevanlinna class function I?A/'(A) for Im A < 0, with 
poles at the eigenvalues of the corresponding Schrodinger operator Lp — Lint ■ The 
relative DN-map is a pseudo-differential operator of order 1: for W|(ri) solutions u 

3/2 1/2 3/2 

the DN-map acts from W2 (P) to W2 (P) and for W2 (il) generalized solutions 
the D-map acts from Wi(r) to L2{T). 
We consider also the boundary problem 

du 

-Au + Vu^Xu, ^|p=Pr, u|g^^^^^p = 0. (12) 

and the operator 

L^ = -Au + Vu, |j|p = 0, Hao.„Ar = 0- (l^) 
with the relative Neumann Green function G^ : 

L^G? = -AG?+VG? = XG?+Six-y), Gf|,^^„^,p = 0, ^|^^^^^^^ = 0. 



The map 

u(x) = J G^{x,j,X)pr{'^)dr =: Q[pp, x e flint , 



(14) 
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gives a solution of the relative Neumann boundary problem (|T2|) . The trace of the 
solution on F 



^ dn 



r 



defines the relative Neumann-to-Dirichlet map which is inverse to the relative 
Dirichlet-to-Neumann map defined above, 

For W2 solutions u the corresponding DN-map acts, on the set of all regular 
spectral points A of the Neumann Schrodinger, from W2^^iT) onto W2^'^{T). For 
W2^^ solutions the ND-map acts acts from L2{T) onto VF2^(r). 
The coefficients of the scattering Ansatz (O can be found, in principle, from the 
infinite linear system which is obtained by substitution of the scattering Ansatz 
into the matching condition (see An important part of the calculation is the 

proof of the formula for the DN-map in terms of the Gp (see or, respectively, 
a similar formula for the ND~map in terms of . Selecting E± as indicated 
above, (O, represent the ND-map of by 2 x 2 operator matrix with elements 
AfV±^± = P±NVrP± 

The similar decomposition of the DN-map of the Schrodinger operator on Vlint 

- ( Zz Z:: ) <«) 

was used in [43j in the course of construction of a convenient representation for 
the scattering matrix on the open spectral bands. We set, in agreement with the 
above notations in (I4I5I6P : 



Hereafter we use the standard bra/ket notations, e >< e' : u — > e < e', u >, with 
the bar on the first factor of the dot-product in i? = L2(r). The exponents of 
oscillating and decreasing modes on the first spectral band spanned by the vectors 
e± G E± are represented as: 

X±e+ = e±'^+-e+ , C-e- = e'^-^e. . 

The matrices 5'™'"' and s™^'", which are defined by the matching of the scat- 
tering Ansatz to the solution of the homogeneous equation on rii„f, constitute 
respectively the scattering matrix - the square table of amplitudes in front of the 
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oscillating modes in open channels {I — 1): 



M 



7n,n—l 

and the table of amplitudes in front of the evanescent modes 



m,n— 1 1, r>2 

The scattering matrix of the junction is represented (see [43 1 and Theorem 12.11 
below) in terms of the matrix elements VAf, NT) combined in aggregates 

M = VM++ - ^ ^ (17) 

N = MV++-NV+^K^j-^jl^—^NV_+ , (18) 

The width 5 of the leads can serve as a small parameter in the course of calculation 
of the scattering matrix. Thin networks, with small S, are characterized by large 
distance between the neighboring spectral thresholds: 

TT^l + iy n^P _ {21 + 1)^2 

One can prove following f43j that, for a "thin junction" , the denominator VM + 

K- is invertible on a major part of a properly selected auxiliary spectral interval 
A, where the DN-map is represented as a sum of a rational function and a regular 
correcting term: 

VAfr = V ^" ^" + /C^i =: VAf^ + /C^, (19) 
A.GA ^^-^ 

The zeros of the denominator VM + K- on A have an important operator- 
theoretic meaning: they are eigenvalues of the Intermediate Hamiltonian. Hereafter 
we consider the rational approximation (|19p and the corresponding rational ap- 
proximation of VM = P^VAfP-: 

VM— = VM^_ + , (20) 

with a regular "error" /C^_ on a complex neighborhood G'(A) of A. 

We call the junction ft thin in closed channels, either in W-^^T) or in W2^'^(r), 
if, respectively, 

sup\\ KZ^IC^_\\wi(r)<'^ , or a^^> \\ KI^IC^_ \\^s,2^^^< 1 (21) 
This implies the following statement (see |35]): 
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Lemma 2.1. // the junction is thin on closed channels, then the denominator of 
( [_?7| ) is invertible 

on a corresponding "major part of the essential spectral interval" - the complement 
of the set of zeros Za G A of the determinant of the finite- dimensional matrix- 
function : 



Za = a : det 



{JC^_+K-) ' VAf^_{\)\ =0} 



1 /2 

A similar statement holds for the above denominator as an operator from W2 (F) 
to W^^^iT). 

Theorem 2.1. The substitution of the scattering Ansatz into the matching con- 
ditions on r gives the following formulae for the scattering matrix on A 

S=[iK++Mr^[iK+-M] , (22) 

S = [JViK+ + [NiK+ - 1] . (23) 

Proof. The scattering Ansatz generated by the entrance vector e € i?+ is con- 
stituted by the incoming wave e'^^+^e, the transmitted/reflected wave e~^^+^Se 
and the evanescent wave e~^^^se: 

The boundary data of the Scattering Ansatz at the bottom sections T should 
match on T the boundary data of the solution of the homogeneous Schrodinger 
equation inside 0™*: 



dip 



= 



= -06(0) = e + S'e + se , 

on 



= Ve(0) = iK+e~iK+Se~K^se . (24) 

r 

Using the matrix representations (fTHl [TS]) for VM , NT), we obtain from (|24|) two 
equivalent linear systems which describe matching conditions on T 

iK+{l - S)e = VU++{1 + S)e + VJ\f+_se , 

-K^se = VN-+{1 + S)e + VU^^se , (25) 

and 

{I + S)e = MV++iK+{I - S)e-MV+^K^se, 
[I +NV^^K^]se ^ NV^+iK+{I - S)eNV^^se. (26) 

Eliminating the component se from them and using the former notations tVJ, M 
we obtain the announced representation for the scattering matrix (|22l23p . 
The end of the proof 

Consider the operator C defined by the above Schrodinger differential expression 
on the junction — flint U lo , with zero Dirichlet condition on the boundary dfl. 
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It is essentially self-adjoint on the domain of smooth functions u, subject to the 
Meixner restriction u e W2^(ri). Assume that the entrance space E — i2(r) on 
the cross-sections F is decomposed as i?_ . We use the former notations P± 
for the orthogonal projections in E onto E±. Consider the Glazman splitting £r 
obtained from £ by imposing an additional partial zero boundary condition on the 
bottom sections F of the leads: 

P+u\v = ' (27) 

complemented by the standard smooth matching condition on F in closed channels. 
The operator C is split by this boundary conditions into an orthogonal sum of two 
operators: 

C — >La®Ia = Ca ■ 

Here I a ~ ^'^j^'^J^^^s in ^^(0, oo) x —: H+, with zero boundary condition 

at the origin u(0) = 0, and La is defined in the orthogonal sum of the channel 
space L^{Q, oo) x E_ 7i_ of the closed channels and L2(^^mt) on W2 - smooth 
functions, subject to the Meixner condition and the matching condition on F in 
closed channels: 

La:D{K) — > L2{n,nt)(sn- . 

Theorem 2.2. The operators La, Ia o-i^e essentially self-adjoint. The absolutely 
continuous components of spectra of the corresponding self-adjoint extensions are 

CoCa) ~ [^1, 00), with multiplicity M, 

00 

cTa(iA) = • (28) 

1=2 l>2 

where each branch cr^ has multiplicity M , and the total multiplicity is growing step- 
wise on the thresholds Xi separating the spectral bands A; = [A;, A;+i]. The spectral 
multiplicity of the absolutely-continuous spectrum of La on the spectral bands A; 
is equal to Ml{l + l)/2. The discrete spectrum of La consists of a countable set of 
eigenvalues accumulating at infinity. The singular spectrum of La is empty. 

The relation MM — I observed from comparison of the formulae (|231 has an 
important operator-theoretic meaning. It is derived from the fact that A^, J\f are 
respectively DN and ND-maps of the Intermediate Hamiltonian — the part L^ of 
the Glazman splitting 

C — > Ca = La®Ia- (29) 

defined by the partial boundary condition ^+u|p = — see [IS]. Contrary to the 
standard splitting C — > Lint ® this splitting (^5]) is finite-dimensional — see 
[3]. The poles of M on the first spectral band, below Xmin, are the eigenvalues of 
A- 
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Figure 6. The intermediate Hamiltonian La inherits the closed 
branches of the continuous spectrum of the unperturbed operator. 
The part of the spht operator inherits the first - open - branch 
of the spectrum of the spht operator. The resonance eigenvalues 
of the intermediate Hamiltonian define the resonance conductance 
of the junction 



3. Krein formulae for the intermediate DN-map 
and ND-map, with the compensated singularities 

Expressions Ai, M va the formulae (|17ll8p contain, at least formally, the singular- 
ities at the eigenvalues of the operators Lp , ■ Presence of these singularities in 
the conditions of the last theorem of the previous section looks strange. In fact we 
were able to prove, see [2], that the singularities in the first and second terms of the 
above Krein formula for M., inherited from L^, compensate each other, so that 
only the singularities of the denominators of (jl7|) play a role. Similar statement 
can be proved, see the theorem 13.31 below, for M . But in fact the compensation of 
singularities permits to obtain more convenient representations for M, Af - that 
is for DN and ND maps of the intermediate Hamiltonian L\. These new repre- 
sentations imply also the corresponding exact formulae for the scattering matrix 
of £, and convenient approximate expressions for the scattering matrix as well. 
This approximate expressions can serve a base for construction of a fitted solvable 
model of the junction in form of star-shaped ID quantum graph, and for derivation 
of the boundary condition at the vertex. 

We begin with the discussion of compensation of singularities in the formula (fT7|) 
for the DN map Ai of the Intermediate Hamiltonian. It appeared, that the sin- 
gularities of the first and second term at the eigenvalues of Lint compensate each 

other, so that only the zeros of the denominator VAf arise as singularities 

of VAf^ on A. A one-dimensional version of the statement can be found in [9] 
and a rescription of the classical Krein formula with compensated singularities is 
given in [44j . In this paper we review the compensation singularities in Theorem 
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13. 1[ following [2J for a general thin junction and prove a similar statement, see 
Theorem 13.31 for the Intermediate ND-map. We also obtain, in course of calcu- 
lations, an important "byproduct" : a version of analytic perturbation procedure 
for groups of eigen-pairs. Note that the standard analytic perturbation procedure 
is aimed on calculation of an individual perturbed eigenvalue. Usually the con- 
vergence of the corresponding perturbation series is limited by the condition of 
non- intersection of corresponding terms As(£) 7^ At(e). Our technique, based on 
compensation of singularities, can be used even to study overlapping terms and to 
study the transformation, under small perturbations, of intersections of terms into 
quasi-intersections. This technique is aimed not on the calculation of an individual 
perturbed eigenvalue, but rather on derivation of an approximate algebraic equa- 
tion for the perturbed eigenvalues and calculation of the corresponding residues at 
the poles of the perturbed DN-map. We also calculate, based on (|^^ . the scatter- 
ing matrix of a "relatively thin" junction in the quantum network. We also develop 
similar technique for ND-map. In fact our technique can be modified to calculate 
the scattering matrix for arbitrary junction, see [2]. 

For given temperature T we consider an essential spectral interval centered at the 
scaled Fermi level A : At = [A - 2m* Th,-^, A + 2m* T tr"^]. We assume that the 
temperature is low, so that At is situated inside the auxiliary spectral interval or 
an open spectral set A 



Our prime aim is: to construct on At a convenient local "quasi-one-dimensional" 
representation for the intermediate DN-map and one for the scattering matrix Sof 
the junction with compensated singularities inherited from the Lint, to substitute 
previous formula ([22]) . 

Selecting an appropriate spectral interval (or just a spectral set ) A : At C A, we 
represent the DN-map PA/" of Li^t on the essential spectral interval At as a sum 



of the rational expression constituted by the polar terms with singularities at the 
eigenvalues € A, s = l,2,...iVof the operator L^nt and the analytic operator- 
function JC^ on a complex neighborhood Gat of At. 

We will also use the operators obtained from VMint via framing it by the projec- 
tions P±, for instance: 



^Compare with the popular one-dimensional formula for the scattering matrix in terms of the 
Weyl function, derived in [52] and intensely used by B. Simon and F. Gezstezy in their approach 
to the spectral inverse problem, see [19) . 





(30) 



P+VN^ntP- = P+VN^P- + P+K^P- = VN^_ + K 



'A 
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We introduce also the linear hull — VfLi {fs} - an invariant subspace of Lint, 
dim E'^ — N, corresponding to the spectrum of Lint contained in A and the part 
L'^ := \s(ps){Vs of Lint in it- To calculate the intermediate DN-map A4 in 

terms of the standard DN - map of Lint we have to solve, see p.ip the equation: 



on the essential spectral interval At- It can be solved based on Banach principle 
if can play a role of a large parameter. 
Definition The junction, for which the operator 



exists on At is called hereafter "relatively thin junction", for the selected spectral 
set A and given temperature T . A general network is called thin, if all junctions 
of the network are thin. 

For a relatively thin junction, due to continuity of JC^_,K^ there exist also a 
complex neighborhood Gat oi At, where IC^_ + K- is invertible. 
Hereafter we assume that the junction is thin. The case of an arbitrary junction 
is considered in [2]. 

The above definition of thin junctions (and networks) is based on the following 
motivation. The DN-map of Lint is homogeneous degree —1. It acts from ^^(F) 

1/2 

to W2 (F), see ^73j. If ^int has a small diameter d then, the norm of the correcting 
term K. is estimated, generically, on the complement of the spectrum, as Const 1/d. 
The same estimate remains true for P_/C'^P_ ;= ]C^_. The exponent on the 

essential spectral band A acts from w'^^'^iT) to W2^^^(F) and the norm of its 
inverse is estimated as Const 5. Then the PFj^^" norm of KZ^ IC^_ is estimated 
generically as Const S/d. Hence, in particular, + JC^_ = if_ [/ + KZ^ /C^_] 
is invertible if S/d << 1, see more comments in 43 . Notice, that for arbitrary 
junction an auxiliary Fermi level Af := Ai can be selected, see 2 , such that the 
condition ([5^ is fulfilled. Now we proceed assuming that is fulfilled. 
Consider the part L^^j of Lint in the subspace E^ = © J2x sa Vs- 



(31) 



-1 



(32) 




A. (^s) {'fs ■■ ^ E^, dimS^ = N. 



Assume 




r 



are linearly independent. Then dim Va eA = dim 



\/ 

VA^eA dn 



N. Denote by T the map 



r 




and introduce 
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It is obvious that dim < J^^^ 



I j<^£s_ I jjj^ Later we will utilize a stronger 
I rJs 



T) s 



Assumption 1 The vectors 
any A G A, hence both : dim < 



W.j{X) :=det<^ (:7 



dn 



are linearly independent in for 
= dim E^ ^ N and 
dipt 



J- 



dn 



N 



(A) > 0. 



(33) 



The above condition ([33]) is equivalent to the pair of conditions: 



1. The functions 



1, 2, 3 ... iV are linearly independent. 



2. The operator J is invertible in the linear hull Vs=i 



s=N dips 



for any A G A. 



Hereafter we reduce the problem of compensation of singularities to the spectral 
analysis of the Schrodinger - type equation 

[L^ - g(A)] ^ = AV' 



in E^ with the A-dependent "potential" 

/ 



Q(A) :=r 



-r+ : E 



A 



E"- 



Lemma 3.1. For thin junction, on the essential spectral interval A the derivative 
^ is a positive matrix N x N . 



Proof. Recall that the branch of the square root a/tt^ 5^'^ + Vg — X is defined such 
that ^^^^ < on the conductivity band. The correcting term /C^ is a meromorphic 
operator - function with a negative imaginary part in the upper half-plane 9A > 

and a positive imaginary part in the lower half-plane, hence — :tt— < on the 
conductivity band Ai. This implies: 



dQ 



r- 



dK^ die 



A 



d\ 



dX 



r+ > 0, for A e At e Ai 



The end of the proof 

Based on some cumbersome calculation, we are able to derive, see [2], that all 
singularities in the Krein formula, inherited from the eigenvalues Ag of the unper- 
turbed operator Lint, are compensated. 

Theorem 3.1 (Compensation of Singularities M.). The Krein formula |j7| j for the 

intermediate DN - map, can be re-written, for a thin junction, on the spectral 
interval A, as: 



M - VU^ = Mreg + J'r+) 



XI^ -L^ + QiX) 
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MA) + ^r+)^^^^^^^^(r^-, (34) 

where — /C^_ JC^^ =: fc(A) is a regular part of M. on At- The 

representation ^34^ remains valid on a complex neighborhood Ga of the spectral 
interval A. 

Remark 1 The announced rescription ((34|) of the Krcin formula (|17|) for the DN- 
map of the intermediate Hamiltonian, has on the essential spectral interval only 
non-compensated singularities, at the eigenvalues of the intermediate Hamiltonian, 
calculated as zeros A J of the denominator A/^ — L'^ + Q{X) := d(A): 

d(A?) = 0. 

These singularities coincide with the eigenvalues of the intermediate Hamiltonian. 
We call the above formula for VJ\f^ the modified Krein formula. Inserting 
into the above formula (j22p gives a convenient representation for the scattering 
matrix of the relatively thin junction, which permits, in particular, to calculate 
sharp resonances, situated near the continuous spectrum, based on eigenvalues of 
the intermediate operator. 

In the case of one-dimensional zeros of the denominator d(Aj)(i^^) = Q.J'^ 
T+i/^ 7^ 0, the corresponding residues are calculated as projections onto the sub- 
spaces 

£f = ^(A?)T+ vf. 

For multidimensional zeros of the denominator, d{\^)N^ = 0, dim N'^ > 1 the 
residues are projections onto the images of the corresponding null-spaces = 

The above expression is analytic in JIa on the complement of the set of zeros 
of the denominator d(A). This means that the eigenvalues of the Intermediate 



zeros of the denominator. 



Hamiltonian are selected from the set. Assume that A^, s — 1, 2, 3, ... are simple 



Theorem 3.2. // the Wronskian \3SI\) does not vanish at the algebraically simple 
(first order ) zero A^ of the denominator, Wj(A^) 7^ 0, then the zero is an eigen- 
value of the Intermediate Hamiltonian, with the same spectral multiplicity. 

Proof. It is sufficient to prove, that the zero is a first order pole of the intermediate 
DN-map, with a finite-dimensional residue having the same dimension as the zero 
of the denominator. Consider the equation 

d{\)u = -L^ + g(A)] u = /, XeE^. (35) 

Assume that d(A^)eJ — 0, and denote by an orthogonal projection onto the 
multiple eigen-space V^, ef of the operator L"^ — Q{Xf) —: Lq, and by R'^ the 
corresponding resolvent: 

- QiX'^) - Xl^r' ^ Rf{X). 
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Then, for A close to we can substitute the "potential" Q in the above equation 
by the Taylor expansion: 

g(A) = QiX?) + (A - X?)§ixf) + ^^^^(A? 

and represent the resolvent near the pole A^ as 

Pi 



R 



Q 



Q 1 



(36) 



(A? - A) 

where i?^ is the part of the resolvent in the complementary invariant subspace 

= (/— Pi)E. Replacing Q(A) by the corresponding first order Taylor formula, 
we rewrite the equation d(A)ii = / by the equation 



^ + i?«(A? 



(37) 



To calculate the residue of the solution u at the pole Af we multiply ([57)1 by ^ 
and by the spectral projection of Lq at A^ and take into account that the 

resolvent of Lq has a simple pole and neglect the terms vanishing at A^, in 
particular all terms arising from the above Taylor expansion beginning from the 
second (A - X^f^{X^). Due to positivity of |2 the operator / + P^^P? is 
invertible. Then, due to (l36t 



P?§u=- 
dX I 



has the polar part at A^ 



A' 



Q 



X 



-/ + 



(38) 



The operator in the square bracket is positive and hence has the spectral form : 

-1 



^1 



I 



pQ pQ 
1 dA 1 



Then the polar term of d^^ = - [Lq(A)] ^ at A? is {X-Xf)-^J2r 
the pole part of the intermediate DN-map at A^ is 



ar Vr) {Vr 



and 



A -A? 



Due to the orthogonality of z/j and non-degeneracy of the Wronskian W^, the 
vectors J^T^ are linearly independent, hence A^ is a simple pole of the inter- 
mediate DN- map, with the spectral multiplicity dim \/^ e^. 

The end of the proof. 
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The scattering matrix of the original problem on the essential spectral interval 
can be obtained via replacement in (|22p the intermediate DN-map by the expres- 
sion (|34p with compensated singularities. This substitution is possible for thin 
junctions, when the exponent in closed channels can play a role of a large 
parameter, compared with the error JC^_ of the rational approximation T^M^ of 
VM. 

This condition may be not satisfied, for given quantum network, at the scaled Fermi 
level A. In that case another representation of the scattering matrix pS)) can help. 
We consider now the problem of compensation singularities for the intermediate 
ND-map N on the essential spectral interval. 

Denote by ips,^^ the eigen-pairs of the operator , see ([T3l) . Select the eigen- 
values from the spectral interval A2, to be defined later, and introduce E^^ := 
Va"gA2 ''^s ^^'^ VA«eA2 ^s|r' ^^'^ consider the map 

t: J2 ^s\r){'^s: Ea,^ E^,s = l,2,...N (39) 

Af eA2 

Assumption 2, see the paragraphs 1 and 2 below: 

l.We assume that the families {V'sjfLi i {V'slfLi |r '^''"^ linearly independent, thus 
are bases in their linear hulls E^^, E^ ' "^^"^ -^A2 ~ '^'"^ ^r' ~ 
Represent the compact in L2{T) relative ND-map NV^^^ —: AfT)^ as 

^ ^ A«eA2 

and consider the corresponding matrices 



and 



-+ 

in the basis of _E = L2(T). Hereafter we omit the upper index A2 on matrix 

elements JC± ±. To represent AfD =: J\f in the form with already compensated 
singularities inherited from the resolvent of , we have to solve the equation 

(I + TVI?— A:_)u = 7VX'_+/ (40) 

Our second basic assumption is the following: 

2. We assume, that the width S of the leads, the essential spectral interval At ='■ 
A C A2 and the rational approximation IC are selected such that 

I + iC^l_ (41) 

is invertible for A G A. 
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For low temperature (that is for a relatively small essential spectral interval At) 
this condition is equivalent to the corresponding condition imposed just at the 
scaled Fermi level A. It is satisfied, if 

sup;,^^ ||<1. (42) 

We will not give here a formal condition which guarantees 2, but just notice that 
due to compactness of the resolvent on for any At there exist A2 D At and 
the corresponding number A^As —■ ^ such that the error of the finite rational 
approximation of the resolvent is small: 

G^ix, s, A) = G^{x, s, /i) + (A - ^I)Q2 '' = 

G (x,.,^) + (A-M)l^ +(A-M) ^1.^ (A,-A)(A,-M)2- 

G^{x, s, fi) + Q2{x, s, A, /i) + K.'^^ix, s, A, ^) =: Q{x, s, A, ^) + K.{x, s, A, ^). (43) 

Here we choose ^ large negative, so that G^{x, s, /i) is a kernel of a small integral 
operator and denote hereafter 

G^ {x, s, /i) + IC"^^ {x, s, A, fi) —: K,{x, s, A, /i), Q2{x, s, A, /i) = Af'D{x, s, A, /i). 

Then for the error P_/CP_ =: /C of the rational approximation Q framed by 

the projections onto the corresponding estimate (|42pis valid. Denote 

{P+ - ]C+^K_{I + ,e__A^_)-i} r+ jf+, 

fK-{I + K.—K-)-^f+ V{\) 
L'^" - XI'^' +V{X) =:i'^^(A), 
r{p+ -i^_(/ + ^__i^_)-i^_+} =:fj+. (44) 

Theorem 3.3. (Compensation of Singularities 

/ 

L'MA) 



= AA.eg - r+ TAr7TT^^+ (45) 



Proof We treat the equation (|40[) with use of the Banach principle under assump- 
tion g2l): 

-ftT-M + X_ (/ + K—K-)-^NV'^^_K^u = K^{I + IC—K^)-^MV^lf. (46) 

Recall that A/"!? = AfV^^ + iC and notice that MT>^^ is connected with the part 
L^^ of in i;'^^ as 

Then, denoting 

7 - 
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we rewrite the above equation (j46p as an equation for v and obtain the solution of 
it in terms of the inverse of the matrix L^'^ - XI^^ + T K_ {I + JC^'^_K~_+)-'^f 
L^^(A): 



Substituting the result into (|46)) 



IA2 _ X1A2 
I 



Then TV/ = 



^A2 _ xi^ 

I 



-T + JC- 



P+f- 



p+r 



TP+f + JC++f-P- 

t 



1A2 _ xi'^^ 

K^{I + IC—K^)-^ 

Leading terms inside parentheses give: 



T 



_ ;^/A; 

I 

L^2 _ xi^ 



P+f- 



-TP- + IC+- 



I 



- 



T + JC 



P+f- 



1A2 _ 



-V- 



T + JC- 



I 



P+f 



L^^-XI^'^ 7^A2(A) lA2_xI^'^ 
I 



K-{i +^c—K-)-^f+- 



L^^{xy 

Taking into account the leading term of the first addendum, we obtain: 



(47) 



l^A2 _ XI^^ 



IA2 _ xi^^^ L^^iX) 



Lower order terms containing 



^ in parentheses give: 



(48) 



I 



I 



lA2_xi^^ L^^(X) 



V 



TP+f^ 



rTP+f, 



LA2(A)--' (49) 
and the adjoint expression. The term which contains only the main singularity 
[LA^(A)]-i is 

L^2{X) 

The terms which do not contain singularities result in: 



:J\f., 



reg 



(50) 
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Note that the operator K^{I + JC K )^^ is self-adjoint on the first spectral 

band. Then, collecting all terms we obtain the announced result. 
The end of the proof 

Based on the theorems 13 . 1 [ 13 .3] we can calculate the scattering matrix either in the 
form (j22p or in the form (j23p . One of these formulae can be more convenient on the 
essential spectral interval, than another, depending on localization of singularities 
oiVN^ andNV^. Luckily, due to to VN^ NV^ = /+, the singularities of the fac- 
tors do not overlap, hence for any point Aq G A one can select an interval centered 
at Ao where at least one of the factors VM^ or MV^ can be substituted by the 
corresponding approximate expression based on the bi-linear formulae suggested 
in [60]. 

4. Approximate Scattering Matrix and the boundary condition at 
the vertex of the quantum graph 

The standard method of calculation of the scattering matrix requires solving of 
an infinite algebraic system anyway, though practically admits a certain simplifi- 
cation in closed channels, see [45] q where it was done for waveguides with simple 
geometry. The approach based on the intermediate DN-map gives a finite linear 
system for the Scattering matrix derived from matching of the component VP^- of 

the scattering Ansatz on the first (open) channel in wires with p = \J \ — Vs — 
on the first spectral band Ai: 

:= e'^'^e^ -f e"'"^ S{p)e^ (51) 

to the limit values on the spectrum, 3A 0, of the solution of an intermediate 
boundary problem with the boundary data on T defined by the scattering Ansatz 
^f. The boundary data ^'j^,, are connected by the intermediate DN-map, 

hence the required matching gives a finite linear system for S: 

ip [e^ - S{p)e^] = MiX) [e^ + S{p)e^] . (52) 

Solving this equation we obtain the formula for the scattering matrix of the op- 
erator C on the first spectral band Ai in terms of Ca by the formula, see (P^ . 
Due to Theorem 13. 21 one can substitute, for a thin junction, the intermediate DN- 
map A^(A) in ([^^ by an approximate expression for A4 — Ai^ + JC^, where 
A4'^ —: A4 approx is a rational approximation VM^ of M on the essential spec- 
tral interval Ay, containing only polar terms with poles on an auxiliary spectral 
interval A, and K,^ is a regular part of on A. 

Theorem 4.1. The resulting approximate expression for the scattering matrix 

S « [ipP+ + M^]-^[ipP+ - M^] Sapprox, (53) 

with p = \/X ~ Vs — , can be used as a first step for the calculation of the 

exact scattering matrix via an analytic perturbation procedure. 



^ The author is grateful to V. Katsnelson for important comments in that connection. 
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Proof Indeed, due to the above theorem 13.21 the error Ai — -Mapprox = A^^, with 
/C^ containing the regular term Mreg too, is real and estimated by 0{6 dj^^). Then, 
due to 13.11 we can represent the exact scattering matrix in form of a product: 

5 = (/ + [ipP+ + M^]-i/C^)~' Sapprox {I - [ipP+ - M^r^K^) . (54) 

Here , K.^ are hermitian on Ai, hence || [ipP+ ± M.^] ^ \\« <5, for thin 
junctions. Hence the analytic perturbation procedure of the calculation the left 
and right factors of the expression in (j54p is geometrically convergent due to 
5 0{5 (i^j-^f) << 1. Thus the exact scattering matrix can be obtained from S approx 
by an analytic perturbation procedure. 
The end of the proof 

One can construct various approximations for the scattering matrix based on 
S approx, replacing M by various approximate expressions, with controllable er- 
rors, see for instance (|58l59p . 

4.1. Simple resonance eigenvalue of the Intermediate Hamiltonian 

The simplest approximate formula for the scattering matrix can be obtained in the 
case when there exist a single simple eigenvalue of the Intermediate Hamiltonian 
on the auxiliary spectral interval A. Indeed, substituting the intermediate DN- 
map A4 — K, + ^Zv^ t)y the corresponding polar approximation generated by 

the resonance eigenvalue Xi of the Intermediate Hamiltonian and the boundary 
current of the corresponding eigenfunction 



2pA_ p d^i 



dn 



dr 



we are able to obtain, due to preceding theorem (|4.1|) . the scattering matrix of a 
thin junction via an analytic perturbation procedure based on the jump-start 



S. 



jump— start 



iK+ + k{X) + a\ ^ 



A- A^ 



iK+ ~ fc(A) - al ^ 



A - A^ 



(55) 



In the case when the first spectral band Ai is the conductivity band , the above 
approximate expression for the scattering matrix can be represented, with P^ = 
P+ e Pi" and pi = VA - TT^S-^ - Vs and A^ « A as: 

^ ipi - k{\) ~ alj:^ ^ 

S jump-start = A + ' , , ^ , x , 9 T^^l ' 
ipi + fc(A) + Ol{j^ 

It corresponds to the one-dimensional solvable model of the junction, obtained 
via attachment an appropriate inner structure to the vertex, see Fig. [T] This ap- 
proximate expression for the scattering matrix can be obtained via imposing on 
the Scattering Ansatx a A-dependent boundary condition at the vertex, see a dis- 
cussion in [43 . Unfortunately this boundary condition does not correspond to a 
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Figure 7. 1-d model of T-junction 



self-adjoint operator, so that it can't be interpreted in terms of Quantum Mechan- 
ics, the same as prominent Wigner boundary condition, see [75] . We are also able to 
represent the scattering matrix with use of a single Blaschke-factor: SapproxiW — 

p{' = p,^ + e^{\)p,^. (57) 

Notice that the scalar Blaschke-factor 0^ is close to —1 on the essential spectral 
interval 

At : {A : |A - < 2m* kT < p~^(A^)} , 

for low temperature T, and it is close to 1 on the complement. For thin junction and 
low temperature the boundary condition can be reduced to Datta-type boundary 
condition. [13j. see below, formula (|64|) represented in terms of boundary currents 
of the resonance eigenfunction of the Intermediate Hamiltonian. 
Once we already developed the compensation procedure based on the represen- 
tation of the Intermediate DN-map in terms of classical DN-map, we can do one 
more step, expressing the approximate scattering matrix (j55|) in spectral terms of 
the unperturbed operator L^nt on the vertex domain Qint , under assumption that 
it has a single resonance eigenvalue Ai E Ay C A, 

Then, for thin junction, the intermediate Hamiltonian also has a simple eigenvalue 
near to Ai. We assume that the major part of the correcting term K.^ in the 
corresponding rational approximation of the DN-map of Lint is defined by a finite 
sum of polar terms and a regular term Mreg — — ^^^a ^j^- /C^^ , see 

m 

s=M T n ^ 

s=2 ^ 



ipjX - Af ) - aj 
ip(X - A^) + a? 
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^ s=2 * 

T^iPi){^, J = P+-IC+-KZ' [l + IC—KZ'r' P-. 

Hereafter we neglect the contribution from higher terms of the geometrically con- 
vergent series 



P-; 
dips 



s=2 



A — A, dn dn 



on ^ on on On 

Then, with only terms containing KZ^ taken into account , we obtain for 
based on Theorem 13.11 an approximate expression for 



lCt-KZ^K,^~ + 



r A-Ai + (^ 

^ ^ an 



— on 



(58) 

For "very thin" junction one can neglect even first order terms containing K_^, 
everywhere, except the expression staying in the denominator, and obtain from 
()58p a simpler approximate formula: 

I 



A-Ai + 



/ dipi 
\ dn 



— on 



[P. 



di£l 

dn 



2 -'I 



where 



A- Ai' 

Q _ 



Q 



P? = ef){e?,ef^a^^P,^ 



ai 



Pa 



dipi 



dn 



I, A? = Ai 



K 



r 



(59) 



P dn Y dr 

This implies an approximate formula for the scattering matrix on the major part 
of the essential spectral interval At El, for low temperature, once the junction is 

thin on the open channel 

p2(A) = |A - (tt^ 5-^ + Vs)\ »|1 II . (60) 



Roughly speaking, on a complement of a certain small neighborhood of the zero of the 
denominator. 
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Notice that arising of the shape of the resonance eigenfunction of the unperturbed 
operator Lint in the corresponding approximate jump-start formula (j6ip for the 
scattering matrix corresponds to folklore observation of physicists, that the eigen- 
functions react to perturbation slower that the eigenvalues. To derive the jump- 
start approximation for the scatterin matrix, denote \/X — tt^ S^^ — Vs p. Then 
we obtain on At, similarly to Theorem 14. II 

ip — k — ^ "^q 

S{X) W H a^^^ ^furap-starf (61) 

In the case when w A we can replace on Aj- the Blaschke factor in front of Pi 
by —1, which implies on At, for sufficiently low temperature: 

Sfump-startW ~ Pi ^ Pi = Soatta (62) 

When modeling the quantum network by a one-dimensional graph, one can at- 
tempt to define a boundary condition at the vertex which implies the scattering 
matrix ([5^ . Indeed, forming the component '5+ of the scattering Ansatz based 
on (|62p. we see the the boundary values of the Ansatz 

satisfy the following boundary condition, similar to one suggested in [13] : 

*+(0) V = 2Pf I/, 4- *+(0)i' ^ 2Piiy. 
ax 

In terms of the components \E'"(0)i^, ^\l/"(0)i^ of the boundary values of the 
Ansatz on the bottom sections r„ of the wires and with use of the components of 
the boundary currents tfi = {4'i}n=i 



P+^ 
dn 



^ ri (63) 

we obtain: (?/^, ^'+(0);^) = 0, £*+(0)t/ || V^, or 



N 
n 

s *+(o)^ „ £ n(o)'^ _ _± *;(o)^ _ _± n(o)^ 



Vi Vi V'? Vi 



N 



(64) 



Example 4: asymmetric T-junction Consider a two-dimensional quantum network 
f2 constructed as a simplest asymmetric T-junction of three straight semi-infinite 
quantum wires width tt/2 attached as shown in Fig. [5] to the quantum well - 
the square Qint :0<a:<7r, 0<y<7ron x,?/-plane. The role of the one- 
electron Hamiltonian on is played the Dirichlet Laplacian. The spectrum of the 



^but not yet equivalent, see an extended discussion below, in next subsection. 
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Figure 8. Simplest asymmetric T-junction 



corresponding unperturbed Hamiltonian on the quantum well is discrete and the 
eigen-pairs e.g. 

2 

Ai = 2 : ipi^i = — sin a; X siny; 



2 2 

A2 = A3 = 5; 2 = — sina; x sin2y and ip2 1 = — sin 2a; x siny; 

TT ' TT 

P5 = Pl,2 + -P2,l. Pl,2 = ^1,2) (¥'1.2 : P2,l = ¥'2,1) (¥'2,1; 



A4 = 8; <^2,2 = - sin 2a; x sin2y; Pg = P2,2 = ¥'2,2) ('^2,2; 

TT 



2 2 
A5 = Ae = 10; (pi 3 = — sina; x sinSy and ifz\ = — sin 3a: x siny; . . 

' TT ' TT 



2 2 
Ae = A7 = 13; 1^2,3 = — sin 2a; x sin 3?/ and 1^3,2 = — sin 3a; x sin 2y; . 



(65) 

are obtained via separation of variables. We choose the basic spectral interval 
A = [4,6], so that there is only one multiple eigenvalue A2 = A3 = 5 of Z/j„t on 
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that interval, and use the approximate formula (|59p for A4. The normal boundary 
current Ji^2 of the resonance eigenfunctions ipi^2 and 1^92,1 is calculated as 



Jl.2 = 



/ 


dtpi,2 








dipi.2 




dn 




9<Pl,2 


V 


dn 


/ 






dn 




dip2,i 




dn 




9<P2.1 


v 


dn 



\ 



J 



( sin2j/ \ 



— smx 



svay 



- - sin 2x 



\ -7 sin?/ / 



(66) 



The spectrum of the unperturbed Hamiltonian on the wires is absolutely- contin- 
uous and has a band structure, with thresholds {P} , ^ = 1, 2, 3 . . . separating the 
spectral bands. The multiplicity of the continuous spectrum jumps up by three 
units on each threshold. We assume that the first spectral band Ai = [4, 16] is 
the conductivity band. The entrance subspace of the open channel is spanned by 
the inferior cross-section eigenfunctions e^j^ — (4/7r)^/^ sin 2x^, s — 1,2,3 on the 
bottom cross-section T^ = [0 < < 7r/2],s = 1,2,3. The entrance subspace 
of closed channels is the linear hull of the superior cross-section eigenfunctions 
^ (4/7r)i/2sin2;a:^, s = 1,2,3, with / = 2,3,.... We wiU calculate the ap- 
proximate scattering matrix (jump-start) of the thin junction based on the ap- 
proximate formulae for (j581 159p . The approximate eigenvalues of the intermediate 
Hamiltonian are found as zeros of the denominator d(A) represented as a 2 x 2 - 
matrix with respect to the basis <y5i.2, V2,i; 

A- 5 



d(A) = 







dn 
dv2,i 



1,21 ^- 

r - 



,KZ' P- 



dipi.2 I \ 
dn \r> 
OV1.2 \ 



, . (67) 

dn Ir^^- dn IT'' ^ dii IT"'^- ^~ dn 

Taking into account only components of the currents in the second spectral channel 



A-5 

8yi,2 
dn 

9^2,1 



,KZ P- 
KZ^P. 



8y2.i I \ 
dn \r> 
aV2,i I \ 

r' 



K. 



dipi 



dn 



,9^1 



dr 



and introducing the following notations for the integrals 

.7r/4 

2 / sin2a;sin4a;(ia; = 2/3 =: a, 
Jo 

Tl/2 

sin a; sin 4a;(ia: — —4/15 =17, a + ^ — 2/5 —: P, 
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we represent the denominator ((67)) and the inverse [d(A)] ^ as 



d(A) = 



A-5 
A-5 



-a (3 



4x/16- A \ -af3 



A-5 



(A-5) , 



1 



TT 



+ 



1 



4V16 - A / a2 + /32 



7 > < 7 



(68) 



A-5 A-5 + (5Q' 
with 

= TT [a^ + /3^]/4 vl6— A- Here is substituted by the contribution 
A/tt sin4x^) (sin 4a;''- from the second spectral branch in the wires. 
If the scaled Fermi- level is 5, then the corresponding multiple resonance eigenvalue 



of Lint is split into pair of eigenvalues Af^ = 5, A^ « 5 



, and the jump- 



4V16-5 

start approximation of the scattering matrix can be calculated in terms of P_|_ 
projections of the boundary currents of the resonance eigenfunctions <^i,2, </'2,i 



/ 



dip 



1,2 



dn 



= P+Jl.2 = 



+ dn 



dn 
dipi.2 



dn 



\ 



^2 



J 



I —-^ sin2a;j'- Jp^ sin 2a;]'- sin 2j/(iri \ 
sin 2a;^ Jp^ sin 2a;^ sin xdX^-^ 
sin2a;^ /p sin 2a;^ sin 2ydr3 / 



sin 2xo 







16\/2/37r2 
V -4/7^ ] 



—j= sin2a;2 V'l 2- 



+ dn 



2,1 



dn 



= P+J2,l = 



Pa 



Pa 



dV2.i 
dn 

dip2.i 
dn 
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sin 2x5'- /p^ sin2a;j'^ sinj/dFi \ 

— ■^sin2a;^ Jp^ sin 2x^ sin 2x c?r2 

y sin2a;;^ Jp sin 2x^ sin y (ir2 / 

/ 16/37r2\ 



sin 2xn 







— = sin 2x2 ^2 1- 
a/tt 



(69) 



4F_ 



\ -16/37r2 / 

The exponent of the first ( open) channel is represented as (A) = \/X" 
where the projection P+ onto the entrance subspace of open channels plays the 
role of unity /_|- = /i + /2 + I3 in and is represented as 

P+ ~ i/n [sin2x]'")(sin2x]'" + sin2x^)(sin2x^ + sin2x;^)(sin2x^] . 

Then taking into account that the contribution /C++ from the major polar part of 



K. 



Ps 
A-8 



IS 



16 sin2x^) 



(sin 2x^ 



A-8 



: sin 2x^) 



7r2(A-8) 



; sin 2x^ , 



and omitting the vector factors 2/-^/?? sin 2x5*-) and (2/Y^sin 2x;j'- on the right and 
left side of the jump-start scattering matrix, 



ipl+ - 



4J3 

A-8 



S. 



jump— start 




P5 

A-5 



A-5+5« 




ipl+ 




P5 




(70) 



For better approximation of the scattering matrix we should use better approxi- 
mation for if _ , /C^ . 

4.2. Symmetric junction. 

The Datta-type boundary condition ([M]) does not coincide with the original Datta- 
Das Sarma boundary condition (|ll2p . suggested in [T3], because the phenomeno- 
logical Datta-Das Sarma boundary condition was suggested for a T-junction which 
is symmetric with respect to the left-right reflection. The resonance concept of the 
conductance permits to derive, for a symmetric junction, the original Datta- das- 
Sarma condition and interpret the phenomenological parameter p. 
Consider a symmetric junction fl consisting of a square (0,7r) x (0,7r), and the 
quantum wires width it/2, attached in the middle of the sides ri,r2,r3,, see 
(|4]). The role of the one-electron Hamiltonian is played by the Laplacian on f2 
with zero boundary conditions. Similarly to above example we assume that the 
electrons are supplied from the second wire, in the first spectral channel. We 
assume that the scaled Fermi level is A =10 and the conductivity band is 4 < 
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A < 16, the eigenvalues of Lint embedded into the conductivity band are Aq = 
5, Ai = 8, A2 = 10, A3 = 13. The corresponding eigenfunctions of Lint are found in 
previous subsection via separation of variables, see (|65p. The role of the resonance 
eigenfunctions is played by (^1^3, (/93 i, with the eigenvalue A5 = Ae = 10. We also 
use the symmetric and antisymmetric linear combinations of them 

2"^/^[V3i^3 + (^34] =: ips, 2"^/^[(^i,3 - ips^i] =: (pa- 

We denote the corresponding boundary currents as 

dn 



T ^ 

on 



Ja 



and consider the projections of them P^Jsym, P+Jasym onto the entrance space of 
the first (open) channel. We assume that the temperature is low, so that the role 
of an essential spectral interval is played by At = [9, 11]. Then the eigenfunctions 
and eigenvalues of the intermediate Hamiltonian can be found based on Theorem 
13.11 taking into account the approximate calculation of the potential Q(A) of L'^(A) 



ith 



Q(A)« ( '^^■^ )D( ^1'^ 



3,1 



f) — ^ dn Ir - dn iT' ^ dn Ir dn 

\ (P K"^ P 2lEhl\ \ Ip "fa,! I p 9ip3,i 



dn Ir - dn Ir' ^ dn Ir - dn Ir' / 

and thus neglecting Q in the case of thin networks. Hence, in the first order ap- 
proximation, the perturbed eigenvalues of the Intermediate Hamiltonian remain 
the same : A^ = A^ = 10. Due to reflection symmetry of the junction there 
are two eigenfunctions of the Intermediate Hamiltonian which correspond to the 
eigenvalue multiplicity 2 obtained based on the Theorem 13.11 The corresponding 
eigenfunctions and the projections of the normal currents onto are respectively 
symmetric and anti-symmetric: 

Va,s - P+^—\P+- 




dn dn 



OLa Ba, Ips = Q^s , = 7 ="565, 



V2 + 7 




1 



and the orthogonal complement in is spanned by the vector 

2/7 1 -jl= i I 



yi+V \ 1 / V2+^ 



1 



with /3 — —2/7. In the first order approximation /C can be substituted by the 
contribution from the nearest eigenvalue A4 — 8, and hence /C++ = fc = 0. 
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The intermediate DN-map is represented, due to Theorem 13. II bv the formula 

h+ 



A - A? A - 10 ' 

where Pa = Cq) (ea, Ps = Cj) (eg. Then denoting by P^ := P+ - [Pa + Ps] =: 
P^ — Pq, we represent the scattering matrix of the symmetric junction as 



aiPa+aiPs 

W^Q - 

P^ 



Here the role of the scalar Blaschke factor 8i in (l57|) is played by the 2x2 matrix 

^_ zp[X~X^]PQ~[alPa+a',Ps] 
ip[X-X^]PQ + [alPa + a^,Ps]' 
The matrix for low temperature is close to —Pq on the corresponding small 
essential spectral interval At centered at A^. Then the scattering matrix is rep- 
resented as 

S ^ Pq - Pq- 

with the pair of complementary projections Pq and Pq, dim Pq = 1, dim Pq = 2. 
Scattering Ansatz on the model graph 

satisfies the following boundary condition at the vertex x = 0: 

^'(0) = (/ + S)e = 2P^e, ^'{0) = ipP+{I - S)e = 2ipPQe. 

Taking into account that dim Pq = 1, dim Pq — 2, we can re- write the previous 
formulae as boundary conditions imposed on the Ansatz: 

P^^'(O) = 0, i/^(0) is parallel to P^e or 

/ J- l'(n^\ n ^i^^) ^2(0) ^^3(0) 

(e ,7/- (0)) = 0, — ^ = — ^ = — ^. (71) 

Ci £3 63 

where P^ -. e^) {e^, = {ei,e^,e^) = (2 + i3'^)-^/^{l, (3,1). This condition 
coincides with the original Datta-das-Sarma boundary condition, see (|l|2p . Our 
analysis reveals the meaning of the phenomenological parameter (3. 



5. A solvable model of a thin junction 

Generally, a Schrodinger operator with non-constant coefhcients or one in a non- 
standard domain rarely admits spectral analysis in explicit form. For qualitative 
analysis of quantum systems the Schrodinger operator often is substituted by a 
solvable model, constructed by the von Neumann operator extension technique, 
[IT] , see for instance [71 [THIS] and an extended list of references in In particular, 
the substitution of the network by a proper one-dimensional graph, with special 
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boundary conditions at the vertices, looks like a convenient tool for the qualitative 
analysis of the Schrodinger equation on the network. Unfortunately the estimation 
of the error caused by the substitution of the network by the corresponding graph is 
difhcult. Shrinking of a "fattened graph" fl^ to one-dimensional graph was studied 
in numerous papers, see for instance [361 137] . The authors considered a compact 
network ils constructed of the vertex domains Qi„ , with the diameter proportional 
to 6", < a < 1 and several finite leads w™, width S, joining them to each other. 
In [66] they developed, based on [Ml El], a variational technique for description 
of the asymptotic behavior of the discrete spectrum of the Schrodinger operator 
on the quantum network of various grades a of thinness. It appeared that the 
(discrete) spectrum of the Laplacian on the compact shrinking "fattened graph" 
tends to the spectrum of the Laplacian on the corresponding one-dimensional 
graph but with different boundary conditions at vertices depending on the speed 
of shrinking: the Kirchhoff boundary conditions at the nodes, in the case of "small 
protrusion" l/2<a<l, or the homogeneous Dirichlet boundary conditions, in 
case of " large protrusion", < a < 1/2, see pJTl, theorems 1,2,3. 
In this section we consider a thin quantum networks with small protrusion a = I, 
assuming diam S << diam Qm- We will construct a quantitatively consistent solv- 
able model of the quantum network, in the form of a star-graph with a vertex 
supplied with inner space and appropriate vertex Hamiltonian. The scattering 
matrix of the properly fitted model serves as a local approximation - on a cer- 
tain "essential" spectral interval A - of the scattering matrix of the original net- 
work. In contrast to the quoted above results for compact networks, where the 
wave-functions are obtained based on the variational approach, we use Dirichlet- 
to-Neumann map DM^ of an intermediate Hamiltonian La, to derive an explicit 
formula for the scattered waves on the original network, see ([22|l . The scattering 
matrix of the star-graph model is obtained via replacement of in (|^^ by the 
corresponding rational approximation on A, based on Theorem 13.11 This defines 
all parameters of the model in terms of spectral characteristics of La. In course 
of construction and fitting of the solvable star-graph model we also define the 
energy-dependent boundary conditions at the vertex for the Schrodinger equa- 
tion on the graph. In the simplest case when only one resonance eigenvalue Aq 
of the intermediate Hamiltonian is present on A, this condition depends linearly 
on the spectral parameter, and is parametrized by coordinates JR/cq, Q'fco of the 
corresponding resonance. 

Note that in [33] an algorithm for construction of the scattering matrix of the 
quantum graph of the scattering matrices of star-shaped elements is described and, 
in [25] a convenient formula for the scattering matrix of the star-graph in terms 
of the boundary parameters at the vertex is suggested. For extended discussion of 
properties of star graphs see [521 130] ■ 

The star-graph solvable model of the thin junction will be constructed as a finite- 
dimensional perturbation of an orthogonal sum of the non-perturbed Hamiltonian 
Ij^ in the open channels and a finite matrix A acting in the inner space of the vertex. 
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We will choose the parameters of the model such that the model scattering matrix 
coincides with the few-pole approximation of the complete scattering matrix 
S. Then the constructed model will be automatically fitted (i.e. quantitatively 
consistent). We assume that the spectral variable is scaled such that original the 
Schro dinger equation on the wires is just —Au + Vsu — Xu. The scattering Ansatz 
of the model in open channels satisfies on the wires to the same equation as the 
scattering Ansatz on the original network 



(72) 

We assume that the Schrodinger equation on the quantum well Clint — ^o, with 
the same spectral parameter, is represented as 

Lqu := —Alt + Vu ~ Xu, 

with a corresponding effective mass iiq. We assume that the Intermediate relative 

DN-map, VJVu, with respect to F, is calculated and the corresponding 

rational approximation is selected. Then the construction of the vertex part of the 
model will be done with a major change of the original Intermediate Hamiltonian. 
For thin networks an auxiliary spectral interval A is selected inside Ai — [j^, 4-p-], 
and hence does not overlap with the continuous spectrum of the intermediate 
Hamiltonian . Only a finite number N of eigenvalues of the intermediate opera- 
tor are situated on A. Then substitution of on A by the rational approximation 
A4a may cause only a minor and controllable error. Now we will prove that there 
exist a finite-dimensional perturbation of the operator 1/^(BA such that the scatter- 
ing matrix of the perturbed operator coincides with Sa ■ The perturbation will be 
constructed via operator restriction-extension procedure applied to the orthogonal 
sum © A, where Ais an N x N Hermitian matrix : Ea Ea, dim Ea = N. The 
parameters of the model will be properly selected to fit the spectral data of the 
Intermediate Hamiltonian on the original quantum network, within the auxiliary 
spectral interval A. 

Assume that the positive matrix A is defined by its spectral decomposition 

Here > are eigenvalues of A, and P_ = i/^,) (v^ are the corresponding orthog- 
onal spectral projections. The eigenvalues and the boundary parameters (3 of the 
model, see below ((82)) . will be defined later, based on comparison of the scattering 
matrix of the model with the essential scattering matrix 

Sa = [iK+ + VNaT^ [iK+ - VNa] ■ (73) 

Consider restrictions of both Ij^ and A to symmetric operators on the corresponding 
domains. The restriction of j^; = is defined on functions vanishing near a; = 0. 
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Then the adjoint operator is defined on -R+ji ^^'^ the boundary form 

of it is calculated via integration by parts: 

JM V) = V) {u, ;>) = {u'io), vm - (c/(o), v'm, (74) 

where U{0), V{0) G and the derivatives are taken in the outgoing direction on 
r with respect to f2i„. 

Restriction of the matrix A is equivalent to selection of the deficiency subspace 
for the given value i of the spectral parameter. Choose a generating subspace 
N^, \/~A'W~ = such that ^^N, n N, = 0, dim N, = d, set Dq = 

{A — il) [E^ G TV.) and define the restriction of the inner Hamiltonian as A — > 
Aq = ^InA- We develop the extension procedure for general N. and fit it later 
based on spectral data of the intermediate operator, see Theorem 4.1, 4.2. In our 
construction N. C i?^ will play a role of the deficiency subspace at the spectral 
point i, dim N. = d, 2d < N and the dual deficiency subspace is 7V_i = - ^ N. . 
The domain of the restricted operator Ag is not dense in E^ , because A is bounded. 
Nevertheless, since the deficiency subspaces A''_^; do not overlap, the extension 
procedure for the orthogonal sum Iq ® Aq can be developed. We will do it here 
with use of the symplectic formalism, see for instance |S3] • In this case the "formal 
adjoint" operator for A^ is defined on the defect N. + N_. :— Af by the von 

Neumann formula : A^e ± ? e = for e G N^.. Then the extension is constructed, 
see lemmas 3.1-3.4 below, via restriction of the formal adjoint onto a certain plane 
in the defect where the boundary form vanishes (a "Lagrangian plane" ). According 
to the classical von Neumann construction all Lagrangian planes are parametrized 
by isometrics V : N- ^ N- in the form 

In case when the deficiency subspaces do not overlap, the corresponding isometry 
is admissible, and, according to [34] there exist a self-adjoint extension Ay of the 
restricted operator Aq. We construct this extension based on the following 

Lemma 5.1. The lagrangian plane Ty in the defect forms a non-zero angle with 
the domain Dq of the restricted operator Ao . 

Proof Indeed, if Ay is the extension, then on the TV it coincides with the re- 
striction of the formal adjoint, and on the domain Dq it coincides with Aq. Then 
assuming that Ty and Dq overlap, we obtain, for some f^J-Ni, v G Ni 

Applying Ay — il to both parts of this equation, we obtain 

hence /-'- = —2iv — 0. 
End of the proof. 
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It follows from the Lemma that, once the extension is constructed on the La- 
grangian plane, the whole construction of the extended operator can be accom- 
plished in the form of a direct sum of the closure of the restricted operator and 
the extended operator on the Lagrangian plane. 

Note that the operator extension procedure may be developed without assumption 
of non-overlapping, see [34] . In particular, the case dimE^ = 1, which is not 
formally covered by the above procedure, was analyzed in [70j independently of 
[34j . The relevant formulae for the scattering matrix and scattered waves remain 
true and may be verified by the direct calculation. We will use this fact in section 
7 below. 

We will use hereafter notations and some facts concerning the symplectic operator 
extension procedure, see Appenxix and references therein. Choose an ortho-normal 
basis in Ni : {/^} , s = 1,2, ... ,d, as a set of deficiency vectors of the restricted 
operator Aq. Then the vectors = a-Ii ■^'^ form an ortho- normal basis in the dual 
deficiency subspace N-i. Under the above non-overlapping condition one can use 
the formal adjoint operator A^ defined on the defect Ni + N_i = Af: 



d 

u = 

s=l 



Y,[Xsfs+Xsfs] ^ M, (75) 
s=l 

by the von Neumann formula, see [3], 

d 

AqU ^'^[-i Xs fs + i Xs fs\- (76) 



In order to use the symplectic version of the operator-extension techniques we 
introduce in the defect a new basis Wf^ , on which the formal adjoint Aq is correctly 
defined due to the above non-overlapping condition: 

2 A-ir^' ' 2i A~iV'' 

A^Wt = W:, A^W- = -W+. 
It is convenient to represent elements u £ J\f via the new basis as 

d 

u = ^[^tw++f,W-]. (77) 

s=l 

Then, using notations X]f=i is,± we re-write the above von Neumann 

formula as 

A + l-t» A 

The following formula of integration by parts for abstract operators was proved in 

my- 
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Lemma 5.2. Consider the elements u, v from the domain of the (formal) adjoint 
operator Aq : 

with coordinates Q. , : 

d d 

Then, the boundary form of the formal adjoint operator is equal to 

J,{u,v) = {A+u,v) - {u,A+v) = (?:;,f_)„ - (f (79) 

One can see that the coordinates of the elements u,v play the role of the 

boundary values {U'{0), U{0), V'{0), V{0)}. We will call them symplectic coordi- 
nates of the element u,v. The next statement proved in [53] is the core detail of the 
fundamental Krein formula [35j, for generalized resolvents of symmetric operators. 
In our situation, it is used in course of calculation of the scattering matrix. 

Lemma 5.3. The vector-valued function of the spectral parameter 

satisfies the adjoint equation [A^ — XI]u — 0, and the symplectic coordinates G 
N. of it are connected by the formula 

Proof see in Appendix, subsection 9.1 or in [53j . 
Introduce the map 

The matrix- function M = Pj^.APj^. — Pjq. Pni has a negative imaginary part 
in the upper half-plane > and serves an abstract analog of the celebrated 
Weyl-Titchmarsh function. The operator- function A4 exists almost everywhere on 

2 

the real axis A, and has a finite number of simple poles at the eigenvalues of 
A. This function plays an important role in description of spectral properties of 
self-adjoint extensions of symmetric operators, see [551 [^ . 

We construct a solvable model of the quantum network as a self-adjoint extension 
of the orthogonal sum lg(BAg. We consider the orthogonal sum of the corresponding 
adjoint 1^ and the formal adjoint: l^ (BA^ , and calculate the corresponding bound- 
ary form J(U, V) :— J{U, V) -\- J{u, v) on elements {U, u) := U from the orthog- 
onal sum of the corresponding spaces. The self-adjoint extensions of the operator 
Z„ © A^ are obtained, based on restrictions of the adjoint operator = © 
onto Lagrangian planes of the form J(U,V). These planes may be defined by 



(82) 
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the boundary conditions connecting symplectic coordinates U'{0), U{Q), , ^" of 
components of corresponding elements in the deficiency subspaces. For instance, 
one may select a finite-dimensional operator /3 : © iV. ^ ® N. and define 
the Lagrangian plane by the boundary condition 

The extension defined by ([5^ on the Lagrangian plane is continued onto the whole 
space {E^ , ) © E^ by forming the direct sum with the closure of the restricted 
operator Aq, see [M]. This construction gives a self-adjoint extension of IqQAq 
in L2{E^, R_^) (B Ej^, defined by the boundary condition ([82]) . 
The absolutely continuous spectrum of the operator coincides with the spec- 
trum of the exterior part of the model, and hence it coincides with the spectrum 
of the trivial component of the split operator £^ (in the open channels). The 
corresponding eigenfunctions of A^ on the first spectral band Ai D A can be 
found, see [4j, via substitution into the above boundary condition for the column, 
combined of the Scattering Ansatz in the open channels with ([50)1 . and, in the 
outer space, with = y/X~Vs~n'^ S^'^: 

^= A+uFu h (83) 

with /3io = /3(J]^ . It gives the linear equation for the Scattering Matrix: 



Solving this equation we obtain the scattered waves and the scattering matrix: 

Lemma 5.4. The scattering matrix for the constructed extension is an analytic 
function of the spectral parameter X: 

with the denominator of the fraction preceding the numerator. The coordinate 
of the inner component of the scattered wave fi8S\) is defined as 



with p — VA— V^+TT^J 



6. Fitting of the solvable model 

It remains to choose the eigenvalues of A, the subspace N. and the matrix param- 
eter /3, such that the operator-function [/J^q -|- I3,,M.(3,„] acting in E^ coincides 

with the essential DN-map T^J^^ of the intermediate Hamiltonian. Denote by 
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the spectral projection corresponding to the eigenvalue k of A, framed by the 
projections Pi onto the deficiency subspace N. 

= PP^P. 

Then the above expression takes the form : 

[Poo+PoiMPJ = 



Poo + ^<^lPoiQ.rPlO 



^ 4 

E^^/^oiQ./?io- (85) 
^ — ' ^ — A 



We will define the boundary parameters = /S^^^ later, but once they are 

defined, we choose /Jqo such that the first summand in the right side of ([85| coin- 
cides with kM /3q(j -|- ^ a Pg-^Q^f3-^g — —kM- Then the scattering matrix takes the 
form: 

4 

^K^-kM+El, ^Jo^QJ^o 

m = , (86) 

iK^ + kM - E„i ^^/3oiQ./3io 

which coincides with the essential scattering matrix if and only if the corresponding 
Krein function 

" 1 * 

kM-Y.^^P,nQJ.o (87) 

■^-^ a — A 

r=l ^ 

coincides with the essential part VAf^ of the DN-map of the Intermediate Hamil- 
tonian on the essential spectral interval At: 

PAA^ « fc(A) + J2 _l ^" • (88) 

r=l ^ 

Summarizing these results we obtain the following conditional statement for the 
extension constructed based on the boundary condition (|82p in case when N. n 
iV_, = or dim ^ 1: 

Theorem 6.1. The constructed operator is a solvable model of the Quantum 
network on the essential interval A, if and only if the dimension of the space 
coincides with the number N of eigenvalues of the intermediate operator on 

A C [Xmax, Xmin], the eigenvalues of the inner Hamiltonian ^ = _i ^r^r) 
coincide with eigenvalues of the intermediate operator on A, there exists a defi- 
ciency subspace N. of the inner Hamiltonian such that Ni n "jz^-^i — '^"'^ 

N 

operator (3^-^ : Ni —* E^ such that for the ortho-normal basis {e^} _^ of eigenvec- 
tors of A in E^ 

^+ ^ = [1 + a']''>oi^«, z.,., r = 1, 2, . . . , iV. (89) 
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Eliminating the inner variables, we can reduce the model to the Schrodinger equa- 
tion with the constant potential on open channels, and appropriate boundary 
conditions on the bottom sections: 



dU 



dx 



\ ' ^+ dn '\'^+ dn 



A - A 



U 



(90) 



Unfortunately, this straightforward construction does not fulfil basic requirements 
of quantum mechanics, and hence we proceed via construction a self-adjoint oper- 
ator in [2.2(0,00) X E^] ® E^). 

Dr. M. Harmer suggested an important strengthening of the previous conditional 
statement, by proving a general theorem of existence of the subspace iV. and the 
projection P^V; which satisfy the condition of Theorem 4.1. The proof we provide 
below only slightly differs from the original proof in 26J: we added an explicit 
formulae for /J^j^, P/y. in terms of the corresponding Gram matrix. 
Denote by the restriction of the intermediate operator onto the invariant 
subspace E^ = E^ corresponding to the part = {Aj , Aj , . . . A„ } of its spectrum 
on the essential interval A, and consider the linear map 



-1/2 



dr 



(91) 



from E. 



to £'^,dim E-^ 



Theorem 6.2. (M. Harmer) The map $^ defines a one-to one correspondence 
between two d - dimensional subspaces, 2d < N : 



A 



E^ ■.= N^ 



If the subspace 



^ ^ and 
is a generating subspace of 



N^n L 



and 
il\ = 0, 



CP, 



E. 



then there exist a unique pair of the boundary operator (3^-^ : E^ 
subspace N. C P^ , which satisfy the condition of the previous theorem. 



(92) 
and the 



Remark 2 This theorem gives an interpretation of the solvable model described in 
theorem 4.1, in terms of the intermediate Hamiltonian via selection of the inner 
Hamiltonian A as a part of in the invariant subspace corresponding to 
the essential spectral interval A. The subspace C P^ plays the role of the 

p! + tl] 



deficiency subspace N. of the inner Hamiltonian and ^P^ — il 
plays the role of the dual subspace N_^. 

Proof The map $^ is represented by the n x matrix $ of columns 4>^ with 
respect to the orthogonal basis of cse-functions {e^} in P^. The condition ((89l) is 
equivalent to the representation of the operator in form P„, , where /J^ is 
a bounded operator acting from P^ into P^ and P„. is an orthogonal projection 
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in Ej^ onto the deficiency subspace TV.. We will construct both /S^^, P„. from the 
data encoded in $a- 

The non-negative Gram operator <i>A $a spectral representation 

$A = D U. 

The non-negative diagonal matrix D is invertible on the orthogonal complement 
E of the corresponding null-space E^^. We denote the restriction D onto E by 
D. One can assume that the subspace E belongs to some extended space E (B Eg 
which contains E^ , and the operator acts from E(BEg onto i?^ as an isometry. 

+ "1/2 

The operator U D coincides with $a- Hence the operator (f>A is presented as 
a product /3P, with f3 = /3,„ = f/^l)''' : E ^ E^ and P = := C i;^, 
dim Ni — dim E — d and coincides with the dimension of the resonance entrance 
subspace of the intermediate operator. Up to some non-essential isometry we may 
assume that E^ — E^ ,A — L^, N. — . The condition ([92| guarantees that 
TV, U = 0. 
End of the proof 

In the case when only one resonance eigenvalue of the intermediate operator 
sits on the essential spectral band, the obtained model scattering matrix 

4 

iK+ - kM + ^f3„,Q,f3,„ 
S(P) (93) 

iK+ + kM - ^^Po^QoP^o 

is a single-pole approximation of the scattering matrix of the network. The condi- 
tion of the above theorem is obviously fulfilled for the single-pole approximation, 
when 7^0, (i = 1, 7V = 1, and (3„ is a one-dimensional operator mapping the 

one-dimensional subspace N. onto the resonance entrance subspace in E^ spanned 

by -P+^f-- For thin or shrinking networks one can estimate, (see [5] and more de- 
tails in [351 H^) the deviation of the single-pole and/or few-poles approximations 
from the exact scattering matrix on the network, in terms of the ratio d/diam 
We postpone the discussion of the non-stationary scattering matrix for QN to 
forthcoming publications. But we notice here that the local wave operators ( see 
[8]) and the corresponding scattering matrix on the essential spectral band can be 
defined for the pair (£, A^). 

7. A solvable model as a jump-start in the 
analytic perturbation procedure 

Recall that the exact scattering matrix was approximated by the essential or ap- 
proximate scattering matrix. In this section we consider this phenomenon from 



46 



B. Pavlov 



the point of view of complex analysis, for the simplest star-shaped network con- 
structed of a single model quantum well with one semi-infinite wire attached to 
it. 

Consider a thin quantum network constructed of a quantum well flin and a single 
quantum wire of width 6 attached to it, 6 /diainQin << 1. Assume that the Fermi 
level is situated on the first spectral band in the wire, which has multiplicity 1. 
Without loss of generality we may assume that the component of the corresponding 
solvable model in the open channel is presented by the Schrodinger equation with 

2/i" = 2/i" ^I,V^ = V,K^=p = ^\-'{ 
E- 



V and one-dimensional subspace 



— u" — p u, < X < oo. (94) 

Assume that the model Hamiltonian is constructed as suggested in the previous 
section based on the "inner Hamiltonian" A, the differential operator and the 
boundary parameters which are reduced to the coupling constant /J^i := /3. Here- 
after we will use the re-normalized eigenvalues — ^ — V := k'^ >0. Introducing 
that notation into the Krein function, and submitting the boundary parameter to 

N 

the condition (3^^ + cttPoi Is Pio = we obtain the corresponding few-pole 

scattering matrix as a function of the wave-number p, with physically meaningful 
limit behavior at infinity S*^ {p) — > /: 



S^(p) = 



ip-k- 13 Y.. 



i+a: 



ip + k + /3^J2 



l+a 



(95) 



wis 



the resonances - sit 
k := —k , when 



Here = |(e,ejp. Zeros p^(/3) of the Scattering matrix ([95)1 
in the upper half-plane Qp > and approach the points ±fc 
13^0. 

Assume that the resonance eigenvalue = k"^ + S^^ + V is situated close to 
the scaled Fermi-level A and the coupling constant (3 := (3^^ is relatively small, see 
below. Separating the resonance term in the numerator and denominator of (j95p 



ip — 13 



1 



+ 



1 



and multiplying by (ip) one can see that — 4- 



plays the role 



of the small parameter. The resonance kg{(3) originated from the eigenvalue k^ of 
the operator A (more precisely : from the point +kg ) can be obtained as a solution 
p = kg {(3) of the equation 



kr, 



P (l + ao)9o 



{p + kg){ip + k + [3 J^^^o p^l^ 



(96) 
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Another resonance originated from the point —A;,-, corresponds to the same eigen- 
value, and it sits at the symmetric point —k^{P) with respect to the imaginary 
axis. Remaining resonances k^{j3), s 7^ 0, can be found from similar equations. All 
functions k^{(3) are analytic functions of /?, k in small neighborhoods of (0, ±fcs). 
They sit in the upper half-plane symmetrically with respect to the imaginary axis 
fc_^(/3) = —k^{P). The scattering matrix (l95l) is unitary on the real axis k and has 
poles at the complex-conjugate points k_, {(3) in the lower half-plane, and hence it 
is presented by the finite Blaschke product which tends to 1 when |fc| — > 00: 

The outer component of the scattered wave is presented as 

=e-'f"-|-S'(/fc)e'f". (98) 

It fulfills appropriate boundary condition at the place of contact with the model 
quantum dot. The inner component of the scattered wave can be obtained from 
Lemma 3.4. 

We explore the model scattering problem for small values of /3. Though the res- 
onances depend analytically on /3, neither the scattering matrix (|95|97|) nor the 
scattered wave depend analytically of (/3,p) on the product of a small neighbor- 
hood of the origin in (3 - plane and small neighborhoods of zLkg in p-plane. The 
analyticity is lost due to presence of the points ±fcj, where the resonances are cre- 
ated at /? = 0: both kg{P) and kg{(3) approach the same point k^ when (3^0. 
The corresponding "resonance" factor of the scattering matrix 

s ^ lp-kM][P + Kim .qq^ 

''^^^^ [p--kM][p+kMr ^ ' 

is non-analytic on {Qp x 17^ x ), though k^ = SRfc^ +i'^k^ is analytic function 
of P due to But the complementary factor of the scattering matrix 



is analytic on (fi^ x fi^ x ) and can be expanded into the power series over 
(3 , m — 0,1,2, ... . Assume now that the function kg{f3) is known. Then the 
following statement is true: 

Theorem 7.1. There exists a one- dimensional perturbation of the operator 



L 



= 



with a non-trivial inner component, such that the scattering matrix of the pair 
(A^ , Ig) coincides with — S^(p). Then the scattering matrix S of the complementary 
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pair (A^,Ajj) is equal to the complementary factor —S ^{p) : 

The complementary factor is an analytic function of (/3, k) on the product (f2^ x 
f2^, X ) of a small neighborhood of the origin in [3-plane and a small neigh- 
borhood of the pair (fc„, — fc^) in the p-plane. 

Proof is presented in [55] • 

Corollary 1 The exterior component of the scattered wave of the operator 
presented by the Ansatz with taken in the form (|57|) is not analytic with 
respect to the coupling constant := /3 near the origin. The non-analyticity of 
the scattered wave is caused by the presence of the non-analytic factor — in 
the scattering matrix. From theorem 4.1, we interpret this factor as the scattering 

matrix for the pair (^A^, l^^. The complementary factor — S° is analytic with 
respect to the coupling constant f3. It can be interpreted as the scattering matrix 
for the pair ^A*^, . Summarizing our observation we suggest, for our example, 
the following two-steps modification of the analytic perturbation procedure on 
continuous spectrum: 

a. First step is the construction of the solvable model and calculation of the cor- 
responding (non-analytic with respect to the coupling constant f3 at the origin) 
scattering matrix. This is "the jump-start" of the analytic perturbation procedure. 

b. Second step is the calculation of the analytic factor of the scattering matrix of 
the model by the standard analytic perturbation procedure. The analytic factor is 
interpreted as the scattering matrix between the constructed solvable model and 
the perturbed operator . 

The obtained connection between resonances and analytic perturbation series on 
the continuous spectrum recalls the connection between small denominators in 
celestial mechanics and divergence of perturbation series, observed by H. Poincare, 
see [62]. More historical comments about intermediate Hamiltonian and the jump- 
start may be found in [55j , where similar modification of the analytic perturbation 
procedure for the Friedrichs model is suggested. 

Remark 3 Note that recovering exact information on the resonance {(3) and on 
the corresponding residue for the perturbed operator A^ , which we need to develop 
the "jump-start" procedure, may be a tricky problem almost equivalent to the 
original spectral problem. On the other hand, if the approximate resonance factor 
is used instead the the exact factor, then the division of the scattering matrix 

through S'' would not eliminate singularity, hence the complementary factor of 
the scattering matrix would be still non-analytic at the origin and hence could not 
be obtained via analytic perturbation procedure. 
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9. Appendix: symplectic operator extension procedure 

John von Neumann in 1933 has found conditions which guarantee existence of 
a self-adjoint extension of given unbounded symmetric operator, and suggested a 
procedure of construction of the extension, see symplectic version in |53j . For given 
symmetric operator Aq defined on Dq in the Hilbert space H, see [3] and given 
complex value A, 3 A 7^ of the spectral parameter: 

Definition 9.1. Define the deficiency subspaces 

Nx ■.= HQ [Ao - XI] Do, 
N-^:=He [Ao-Xl] Do. 

The dimension of N\, is constant on the whole upper and lower spectral half- 
plane 3A > 0, 3A < respectively. 

Definition 9.2. Introduce the deficiency index {dimNx, dimN^) :— n_) of the 
operator Ao ■ 

J. von Neumann proved that 

Theorem 9.1. The hermitian operator Ao has a self-adjoint extension if and only 
if n+ — n_ 

The idea of construction of the extension is based on the following theorems von 
Neumann, see for instance [3j: 

Theorem 9.2. The domain of the adjoint operator is represented as a direct sum 
of the domain D^^ of the closure and the deficiency subspaces, in particular: 

D^^=Da,,+N, + N^,. 

The deficiency subspaces of the densely- defined operator are the eigen-spaces of the 
adjoint operator: 

Theorem 9.3. If Ao is an Hermitian operator with deficiency indices (n+, n_), 
n_ = and V is an isometry V : Ni ^ N-i. Then the isometry V defines 
a self-adjoint extension Ay of Ao, acting on the domain 

Dav = Dao + {e» + Ve„ a G A^;} 

as a restriction of Aq onto Dav ■ 

Ay : uo + ei-{- Vci — > ^0^0 — isi + iVe^. 
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J. von Neumann reduced the construction of the extension of the symmetric opera- 
tor Aq to an equivalent problem of construction of an extension of the correspond- 
ing isometrical operator - the Caley transform of ^o- It is much more convenient, 
for differential operators, to construct the extensions based on so-called boundary 
form. 

Example 5. Symplectic Extension procedure for the differential operator Consider 
the second order differential operator 

defined on all square integrable functions, u G L2(0,oo), with square- integrable 
derivatives of the first and second order and vanishing near the origin. This opera- 
tor is symmetric and it's adjoint Lq is defined by the same differential expression 
on all square integrable functions with square integrable derivatives of the first 
and second order and no boundary condition at the origin. This operator is not 
symmetric: its boundary form 

J{u, v) = {L+u, v) - {u, L+v) ^ u'(0)5(0) - ii(0)w'(0), u,v £ 

is generally non equal to zero for u,w e D^+. But it vanishes on a "Lagrangian 
plane" V-^ C defined by the boundary condition 

it'(O) — 7u(0), 7 = 7. 

The restriction of the Lq onto the Lagrangian plane is a self-adjoint operator 
in L2(0, 00): it is symmetric, and the inverse of it [L^ — XI)~^, at each complex 
spectral point A, exists and is defined on the whole space £2(6, 00). 
The operator extension procedure used above for the differential operator, can 
be applied to general symmetric operators and serves a convenient alternative for 
construction of solvable models of orthogonal sums of differential operators and 
finite matrices. We call the abstract analog of the extension procedure the sym- 
plectic extension procedure. Let A be a self-adjoint operator in a finite-dimensional 
Hilbert space E, dim E ~ d, and iV,; is a subspace of E, dim N ^ n < d/2, 

which does not overlap with jV^ :— N^i: 

n iv_« = {0} . 

Define the operator Aq as a restriction of A onto Dq :— -j^EQ N . This operator 
is symmetric, and the subspaces N±i play roles of it's deficiency subspaces. The 
operator can Aq can be extended to the self-adjoint operator Ay D Aq via sim- 
plectic extension procedure involving the corresponding boundary form: selecting 

Introduce, following [53|, another basis in the defect N ^ Ni + N-i 

A + il 



2 



A + il 



2i 



A-il 



9s 



Due to A^gs + igs = 0, [A^ - il]^%9s = we have. 
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Following [24] we will use the representation of elements from the domain of the 
adjoint operator by the expansion on the new basis: 

s 

with uo G D{Ai)) and simplectic coordinates ^4- 

We also introduce the boundary vectors of elements from D[A'^) 

4 eA^^, 



Define the formal adjoint operator A^ on the defect N = A^^ + N^i as: 

AqC^ = — «e+, fore-|_ G A^^, A'^e^ = ze_, fore+ G A^-j, 



Then we have: 



A+W+^W-, A+W-^-Wj 



Following [53j , we will use the representation of elements from the domain of the 
adjoint operator by the expansion on the new basis: 

s 

with uq G D{Ao) and simplectic coordinates We also introduce the boundary 
vectors of elements from D{Aq) 

4 :=^?i<?ae7V„ 



u = uo + - -J^C ■■= no + n\ uo G D{Ao) G A^. 

Then the boundary form of Ag is calculated as 

{A+u,v) - {u,A+v) := J{u,v) = (Ct^t) - it, CI) 
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9.1. Operator Extensions: Krein formula 

Theorem: Krein formula Consider a closed symmetric operator Aq in the Hilbert 
space Ti., obtained via restriction of the self-adjoint operator A onto the dense 
domain D{Aq), with finite- dimensional deficiency subspaces A^ip,;, ■— P+, dim 
Ni — dim N^i. Then the resolvent of the self adjoint extension Am defined by the 
boundary conditions 

i+ - Mt (101) 
is represented, at regular points of A]\j, by the formula: 



{Am -\ir' = 



I 



A + il 



-P+M 



I 



P+ 



A-il 



A- XI A-XI ^ i + p.I±^p,M ^A-XI 



Proof. For the convenience of the reader we provide below the sketch of the proof 
of the Krein formula via simplectic operator extension procedure. Solution of the 
homogeneous equation (A~^ — XI)u = / is reduced to finding uo,^± from the 
equation 

(A-A/).,-^e-^e^/, (102) 
Applying to this expression the operator -jEtti due to {A — iI)uo-LNi, we obtain 



A-il 
A -XI 



Then, from the above equation p02p and from the boundary condition poip . we 
derive: 



1 



Mo 



A-il 



and 



I + XA 
A -XI 

A 



A-XI 



f, 



A -XI 
The end of the proof 



f 



Uq 



A + il 

A-xr 



A-il 



-P.M- 



A-il 
I 



I + P- 



A-il 

A-xr 
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